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Abstract. We investigate adiabatic oscillations for S Scuti star models, taking into account a moderate rotation 
velocity (~ 100 kms -1 ). The resulting oscillation frequencies include corrections for rotation up to second order 
in the rotation rate including those of near degeneracy. Effects of either a uniform rotation or a rotation profile 
assuming local angular momentum conservation of the form fi = fi(r) on oscillation frequencies are compared. 
As expected, important differences (around 3 /iHz) are obtained in the g and mixed mode regions. For higher 
frequency p modes, differences range between 1 /Mz and 3 (Mz. Such differences are likely to be detectable with 
future space missions such as COROT, where precisions in frequency around 0.5 jMz will be reachable. 
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1. Introduction 

Intermediate mass stars are characterized by a convec- 
tive core and a radiative envelope. As representative of 
such stars, 5 Scuti stars are located in the lower part of 
the Cepheid instability strip, with spectral types from A2 
to F0. Such pulsating stars seem particularly suitable for 
determining the extent of the convective core and inter- 
nal rotation rates, and thereby to understand better hy- 
drodynamical processes taking place in stellar interiors. 
Particularly, balance between rotationally induced tur- 
bulence and meridional circulation generates mixi ng of 
chem icals and redistribution of angular momentum 
1992), which affects the rotation profile and the evolution 
of the star. It is expected that inter mediate mas s stars do 
not rotate uniformly as a solid body. IZahnl jl992j) proposed 
that, as a result of strong anisotropic turbulence, the stel- 
lar rotation profile is shellular, i.e. the star is divided in 
differentially rotating concentric shells. 

S Scuti stars can be found in two principal evolu- 
tionary stages: on the main sequence, where they burn 
hydrogen in their convective cores, and before the sub- 
giant phase, burning hydrogen in shells. They can be suf- 
ficiently evolved that they present a strong gradient of 
chemical composition, produced by their shrinking con- 
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vective core. These types of structures generate partic- 
ular modes, showing a dual behaviour. These modes are 
known as mixed modes, which behave as pressure modes (p 
modes) at s urface, and as gravity modes ( q modes) toward 
the center l|Christensen-Dalsgaardlll998ll . Indeed, they 
can penetrate sufficiently deep toward the center 
of the star and at the same time show amplitudes 
at the surface large enough to be detected. The ex- 
istence of such modes can thus be very important 
when investigating deep interiors of stars. 

In the last decades, great efforts have been employed 
to develop observational seismology of S Scuti st ars, for 
instance the World-w ide observational camp aigns JBregei 
200rl lHandieil 1200(1 or STEPHI network ^Michel et al. 
2000). However, several observational aspects of the pul- 
sating behaviou r within the instability strip are not fully 
understood ( see iTemnleton et al.lll9*9 7). Due to the com- 
plexity of the oscillation spectra of S Scuti stars, the 
problem of the identification of modes remains unsolved. 
As 6 Scuti stars are commonly fast rotators (100 < 
usini < 200 kms -1 ), additional uncertainties to their al- 
ready complex oscillation spectra arise from the effect 
of rapid rotation. Rotation alters the internal structure 
of a star through a modified hydrostatic balance and, 
likely more important, through mixing cause d by circula- 
tion and /or instabilities induced by rotatio n ||Zahn|[l992; 
iMaeder & MevnetJl200fil iHeger et al.ll2000lb Furthermore 
the simple characteristic pattern of symmetric multi- 
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plets split by rotation is broken. In the framework of a 
perturbation analysis, second-order effects in duce strong 
asymmetries in the splitting of multiplets llSaiol Il98ll: 
IPziemhowski fc Goodd llflfla hereafter DG92) and fre- 
quency shifts_whkhcannot be neglected even for radial 
modes l|Soufi et alJll99fi|) . 



2. Non-degenerate oscillation frequencies 
corrected for rotation effects 

Following DG92, when Coriolis and centrifugal forces are 
considered, the eigenfrequency li and associated eigen- 
functions £ must satisfy the following oscillation equation: 



In DG92, the authors propose a second-order for- 
malism to study the eigenfunctions for the Sun and S 
Scuti stars, with both radial and latitudinal rotation, 
which takes into account the effects of near degeneracy 
up to first order in the rotation rate. For medium-high 
rotators, like 5 Scuti stars, the effects of near degener- 
acy up to second and third order cannot be neglecte d 
(see iDziembowski fc Goodel Il992t iGounil et all 120001) . 
This last point is included in a co mplete third order 
formalism given hv ISoufi et alJ l|l998l) . hereafter SGD98. 
Near degeneracy affects the asymmetry of multiplets, 
making even more difficult the interpretation of ob- 
served oscillation spectra. Based on this approach, 
assuming a uniform rotation for the sake of simplicity, 
theoretical and quan titative studies on oscillations of 
rotat i ng 6 Scuti stars iGounil et a,l.ll2000tlGo ur>il & Talon 
2002;; iDaszvhska-Daszkiewicz et alJ l20 02t IPamvatnvkri 



2QQ3!) ; on solar-li ke sta rs IPzi^mbow^kTfc^oupT 
199 81 iGoupil et alJ 12004) and Cephei stars 
l|Da,szvns ka,-Da,sz kiewicz et alJ 120031) have followed. 
SGD98's work has been recently revisited and applied to 
the s tudy of oscillations of (3 Cephei stars l|Ka,ra,mi et alJ 
1200.4 in press). 

The present work focuses on the magnitude of the 
effect of shellular rotation on adiabatic oscillations of a 
1.8 M© star. In order to avoid possible interferences be- 
tween third order effects of rotation on oscillations and 
those coming from shellular rotation, only second-order 
terms are consider ed. To do so, we have built a numerical 
code JSuarezl2002h taking into account a complete formal- 
ism up to second order in presence of a shellular rotation 
n = ft(r), 

Theoretical oscillation spectra computed for two dif- 
ferent models of S Scuti star are compared: one model has 
been evolved assuming uniform rotation (global conserva- 
tion of the angular momentum), and another model has 
been evolved assuming a shellular rotation (local conserva- 
tion of the angular momentum) . The impact of a shellular 
rotation on oscillation frequencies is discussed in the frame 
of the future space experiment COROT (launch in 2006). 

The paper is organized as follows: Sect. briefly re- 
minds the basis of the second-order perturbation formal- 
ism of non-degenerate oscillation frequencies. Section 
describes near degeneracy theory as implemented in our 
oscillation code. In Sect. 0] the adopted methodology is 
detailed. Section discusses the results and finally con- 
clusions are given in Sect. El 



£o£ - P^ 2 ~ 2p uSlK£ + CG a - 92^ 2 ) £ = 0, (1) 

where Co = lu + mQ and K = 8e 2 x (SGD98, Eq. 22). The 
linear operators £jq and £2 respectively correspond to Lq 
and L 2 defined in SGD98 (Eqs. 23 and 24). The mean 
density of the pseudo rotating model is represented by 
po. Respectively, its perturbation by the non-spherically 
symmetric component of the centrifugal force is given by 
the following expansion in Legendre polynomials 



92 = P22{r) P 2 (cos(9), 



(2) 



where P22M is defined in SGD98 (Eq. 15). With appropri- 
ate boundary conditions, Eq.^ forms an eigenvalue prob- 
lem (appendix EJ. The eigenfunction £ is written as 



(3) 



where £0 and £1 respectively correspond to the non per- 
turbed and first-order perturbed eigenfunctions, defined 



£0 
£1 



2mVt 



to 



yiY e m e r + ziV H 17 



t £+ i e r x Vh^+i + h-i e r x VhY^i 



(4) 



(5) 



where the notations are the same as in SGD98. In these 
expressions, yoi and zq respectively represent the normal- 
ized radial and horizontal co mponents of the ei genfunction 
for the fluid displacement l|Unno et alJ 1989). Similarly, 
the first-order corrections are represented by y\ and z\ 
respectively (see Annendix lA.lJl . 

Hence, the eigenfrequency associated with each mode 
is labeled with the subscripts n, I, m representing, respec- 
tively, the radial order, degree and, azimuthal order of the 
corresponding spherical harmonics. The temporal depen- 
dence is chosen of the form exp(iujt), so that prograde 
modes correspond to m < 0. Following DG92 oscilla- 
tion frequencies are obtained by means of a perturba- 
tive method taking into account up to second order ef- 
fects of rotation. Equilibrium and oscillating quantities 
are expanded with respect to e = fla/uj and /i = il/a, 
where a is the dimensionless mode frequency defined as 
a = uj/iGM/R 3 ) 1 ' 2 ; R and M represent the stellar ra- 
dius and the mass of the model respectively. 

Ignoring in this section the resonant interaction due to 
near degeneracy, the frequency of a given rotationally split 
multiplet of degree i and radial order n can be written as: 



(6) 
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where the subscripts n, £ have been omitted. The w m fre- 
quencies are obtained by means of a perturbative method. 
The zeroth-order contribution, w , represents the mode 
frequency with all effects of rotation ignored except that 
of the horizontally averaged centrifugal force in the equi- 
librium model (later on, we will abusively refer to it as 
the ' non perturbed' frequency). The other two terms, to\ _ m 
and u>2,m represent the first- and second-order corrections, 
respectively. 

A shellular rotation is considered with a rotation pro- 
file defined as 



n(r)=fi[l + »»(r)], 



(7) 



where ^ represents the rotation frequency at the stellar 
surface. In this context, the first- and second-order fre- 
quency corrections are written as: 



wi,m = mQ (C L - 1 - Jo) 

^2 

w 2 , m = — (-Do + m 2 D 1 ) . 

LU 



(8) 
(9) 



The Ledoux constant, Cl, which determines the usual 
equidistant splitting valid in the limit of slow rotation, 
is given by the well-known integral expression: 



1 



R 



Cl = j- I [22/oi^o + zo 2 } Po r 4 dr , (10) 

J Jq 



where Iq is the term of inertia given by: 

Io= [yoi 2 + Az 2 ]p Q r 4 dr. 
Jo 



(11) 



with A = £(£ + 1). In Eq. [HI Jo represents an additional 
contribution in the case of shellular rotation and is given 
by: 

1 f R 

Jo = -r / Vo(r)[yoi 2 + A^o 2 - 2y iz - z 2 ] p r 4 dr . 
h Jo 

(12) 

The second-order coefficients Dq and D\ in Eq. Intake into 
account the non-spherically symmetrical distortion due to 
the centrifugal force. The symmetry of split multiplets is 
broken by the m 2 dependency (see Eq. EJ. For later use, 
we also rewrite these coefficients in Saio's notation: 



D = X x + X 2 



D 1 = Y 1 + Y 2 



(13) 



Comparing Eq.Elwith the formulation given in DG92, the 
following relations are obtained: 

— ) {X 1 + m 2 Y 1 ) =LU J + LU l + ^+^ , (14) 



for X\ and Y\ and, 



wo ( — ) (X 2 + m 2 Y 2 ) = loV 

UJq 



(15) 



for X 2 and Y 2 . The analytic expressions for Xj and Yi are 
given in appendix El The terms u> 2 , u>\ and u> 2 , defined 



by Eqs. 16-19 in DG92, represent the different contribu- 
tions to second-order frequency corrections in absence of 
degeneracy effects and, as in SGD98, include the effect 
of the symmetrical component of the centrifugal force on 
the equilibrium model: uj 2 and usj represent, respectively, 
the poloidal and toroidal components of the frequency, ob- 
tained from the first-order correction to the eigenfunction; 
wJ? corrects for the effects of the centrifugal force on the 
stellar structure. 

3. Near-degenerate oscillation frequencies 

The perturbation method presented in the previous sec- 
tion assumes that the non-perturbed eigenmode can be 
represented with one single spherical harmonic. This is 
no longer valid in the case of near-degenerate frequencies, 
i.e. when two or more frequencies are close to each other 
(u n [ m ~ u> n i m >). Consequently, the perturbation method 
must be modified in order to include corrections for near 
degeneracy as done in the next sections. 

3.1. The oscillation equation in presence of degeneracy 

The formalism is similar for 2 or 3 near-degenerate modes. 
For the sake of simplicity, the formalism is illustrated be- 
low for two modes a and b. The a and b subscripts repre- 
sent (n,£,m)a and (n,£,m)b respectively. For near degen- 
erate mode a and mode b, it is convenient to define: 



W0, a + Wo,f> 



(Wo — W 0j a — Wq, b ■ 



(16) 
(17) 



The eigenfrequency and the eigenfunction of a near- 
degenerate mode are then assumed of the form: 



LU = LUq + U>1 + UJ 2 



(18) 
(19) 



j=a,b 



First and second-order corrections to the eigenfrequency 
in presence of near degeneracy are represented by oj\ and 
uj 2 respectively; £o,j an d £ij are the non perturbed and 
first-order perturbed eigenfunctions respectively defined 
in EqH& Eq.H 

Let us now replace u> and £ by Eq. El and Eq. El 
in the oscillation equation (Eq. Q . Perturbation is then 
performed keeping in mind that Sujq is small, of first or 
second order in the rotation rate i.e. Sujq ~ 0{fl) or Suio ~ 
0(£1 2 ). Projecting onto the non-perturbed eigenfunctions 
£ a an d £o b> the- following system is obtained: 



E 

j=a,b 



M 



(i) 

jk 



SuJr 



3. 



fop j 
8w n 



= 



(20) 



for k = a,b. In this equation we have defined the unit 
matrix 3i and also the matrix 3_i as: 







(21) 
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Fig. 1. Radial component of the centrifugal acceleration 
for a differentially rotating stellar model with Jl (contin- 
uous line), and a uniform rotating one with Jq (dashed 
line), considered in the present study. These quantities 
are displayed for a 1.8 M© model, aged of 1050 Myr with 
a photometric radius of R — 2.27 Rq. 

The first- and second-order matrix M.jk terms are given 
by: 



M (i) = (-th+vxj 

3 \ w l,kj 



-U) 2 + Mj 

where we have defined 



^2 + Mfe 



2 



— (Jij + CJ. 



J 2,jk 



J 2,jk ■ 



(22) 

(23) 

Ai ( 24 ) 
(25) 



Terms with one single subscript are not affected by near 
degeneracy and are defined in Sect. [21 Coupling terms for 
degeneracy are included in wi.jk and fijk; they are defined 
with double-subscript terms and their definitions and de- 
tailed expressions are given in Appendix [O 

3.2. Selection rules 

Corrections for near degeneracy exist only for modes with 
degrees £ and azimuthal orders m fulfilling specific selec- 
tion rules: 

• The analysis of wi, a b (Eq. IC.7II shows that frequency 
corrections for near degeneracy already arise at first or- 
der whenever the modes have same degree £ and same 
azimuthal order m, i.e. £ a — £b and m a — mt. The 
modes then differ only by their radial orders n a ^ rib, 

• Near degeneracy affects second-order frequency correc- 
tions of near-degenerate frequencies (see Appendix [C|l 
whenever £ a — lb or £ a — £b ± 2 and m a = nib, These 



selection rules are also obtained in the case of the com- 
pound third order treatment given by SGD98. Eq. EDI 
is then solved for j = a and k = b for each case allowed 
by the selection rules. 

For two near-degenerate modes, a and b, the proximity 
in frequency is expected to be less or equal to the rotation 
frequency of the stellar model (|w — < ft or fl 2 /ui a ) de- 
pending on the degrees and azimuthal orders of the modes. 
This estimate may vary with the nature of the modes (g 
modes, p modes or mixed modes). For stellar models of 
our interest here, numerical applications reveal that many 
modes have near-degenerate frequencies. 

3.3. First-order near degeneracy 

In this case 5u>q is 0(fl). Only modes a and b with £ a = lb 
modes are affected. Such a situation generally concerns 
only a few modes which are in avoided crossing. We find 
the first-order correction to the frequency as the condition 
for the existence of non-trivial solutions of (Eq. I20|l . 



SuJo 

-U)i+LJi t a+—^- 



-—)-Whab 



where 



Wi, a6 



Wl.ab ba ■ 



0, (26) 



(27) 



In the case of two degenerate modes, this system is equiva- 
lent to Eq. 59-60 in DG92. Conditions for non-trivial solu- 
tions give the first-order frequency corrections in presence 
of near degeneracy 



Ula + Ux t b f7T? 

^1 = -Z ± \/-H-i,ab ■ 



(28) 



with 



2ti,ab=y — — + ~) ( 29 ) 

Note that, in case of negligible near degeneracy, that is, 
when Auj\ ab « {ui\, a — &i,b — Sloq) 2 , the non-degenerate 
u>i t a and uj\,b frequencies are retrieved. The second-order 
near-degenerate frequency correction is again obtained as 
the condition for a non-trivial solution of Eq.20: 



2 


SUJQ 

where we have defined 



6uj 2 \ / 

V a + Ha - U 2 + 7— I ( V b + Ait. - ^2 + 



8Q ) 



"W 2 ,a6 = (Ul,ab + (l>ab)(ui,ab + ^ab) 



w 2 , ab = o. 

(30) 

(31) 



and 



Vb = 



OJl + UJ\h 



2 J 



(32) 



with uji given by Eq.26. Eq. |2H1 provides two first order 
solutions, uif and oj^ . Four possible solutions for cl>2 are 
thus obtained. Although all of them are mathematically 
valid, only two have a physical meaning. Each second- 



order solution 0J2 



and cD is associated to one and 
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Fig. 2. Evolutionary tracks of 1.8 M models representa- 
tive of a typical S Scuti star. The three tracks correspond 
to: a uniform rotating model (dashed line), a shellular ro- 
tating model (continuous line) and finally, a non-rotating 
model (dash-dotted line). The box represents typical ob- 
servational errors for S Scuti stars. 



Fig. 3. Rotation frequency through the star normalized 
to its surface value. The rotation profiles correspond to 
Bl model (curved line, shellular rotation) and Bq model 
(horizontal line, uniform rotation) respectively. Models 
are built with the same surface rotation frequency £lg = 
9.98 /J3z. For more details see Sect. 21 



only one first-order solution Cj 1 and , respectively. We 
then obtain: 



( Vh + Vg Pb + Pa 6uj%\ /— (TT 

" 2 = [—2- + —2~ + ± v 



where 
with 



h(v,n) a ,b 



54i = %,M) 2 a )6 + W 2 , ab 

(v a + fi a ) - (y b + fi b ) 



(33) 
(34) 
(35) 



Table 1. Characteristics of the three 1.8 Mq models con- 
sidered in Fig. |2 From left to right, il s represents the 
surface rotation frequency, T e g the effective temperature, 
X c the central hydrogen fraction and p, the mean den- 
sity. Two evolutionary stages are considered: the zero age 
main sequence (A), and a main sequence age of 1050 Myr 
(B). A solar metallicity, a mixing length parameter of 
oi m l — 1-614 and, an overshooting parameter of dov = 0.2 
have been assumed. 



3.4. Second-order near degeneracy 

In this case, Sujq. is 0(Q 2 ), Q\ — w\ and, Wijfc = in 
Eq. [21 The quantity fij in Eq. ls then u>2j as defined 
in Eq. Only modes with £ a = ib ± 2 are concerned. The 
analytic expression of second-order frequency corrections 
is obtained as the condition for existence of non-trivial 
solutions for 



j=a,b 



We then obtain 



U)2 



Mfe + Pa 



3-1 



0. 



r (2) 



L 2,afc > 



where 



v (2) _ ( Pb - Pa 
JX 2,ab — \ n 



SlOq \ 2 



pab pba 



(36) 



(37) 



(38) 



When near degeneracy is negligible, ^fi a bPba << (p-a — 
fib + Sujq) 2 , and the non-degenerate solutions fi a = u 2a 
and fib = ^2b are respectively retrieved. 



Model 


n. 


R 


log T eff 


x c 


P 




(fMz) 


(Re) 


(K) 




(gem -3 ) 


Al 


15.485 


1.824 


3.870 


0.726 


7.308 


A G 


10.902 


1.800 


3.873 


0.726 


7.555 


Ao 





1.795 


3.876 


0.726 


7.927 


B L 


9.980 


2.271 


3.888 


0.318 


2.839 


B g 


9.757 


2.275 


3.888 


0.315 


2.827 


Bo 





2.242 


3.892 


0.312 


3.939 



Our numerical computations suggest that for rapid 
enough rotations, every pair of modes having A£ — ±2 
with radial orders n and n — 1 respectively, should present 
frequencies near enough to be coupled. For pairs with 
I = 0, 2, coupling is more significant for high radial or- 
der modes. 



4. Evolutionary models 

Equilibrium stellar models have been computed with the 
ev olution code CESAM. Followin g the approach described 
in iKinnenhahn fc Weigertl l|l990l) , equilibrium models are 
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constructed by modifying the stellar structure equations 
so as to include the spherical symmetric contribution of 
the centrifugal acceleration, by means of an effective grav- 
ity. 

g eS = g-Mr), (39) 

where g is the local gravity component at a radial distance 
r 

m(r) 



g = G- 



(40) 



with m(r) the spherical mass at radius r, G the gravita- 
tional constant and, A c {r) corresponds to the centrifugal 
acceleration at radius r : 



A c (r) = -rn 2 (r). 



(41) 



This spherically symmetric contribution of the rotation 
does not change the shape of the hydrostatic equilibrium 
equation. This is not fully representative of the structure 
of a rotating star. 

First, the non-spherical component of the defor- 
mation of the star is not considered. However, its 
effects are included through a perturbation in the 
oscillation equations (see previous sections). Second, 
rotationally-induced mixing and transport of angular mo- 
mentum can sign i ficantly modify the rotation pr ofile 
jTalon et al.lll99?t iMaeder fc Mevned l200fi iHeger et alJ 
200(j iDenissenkov fc ToudbnOft IPalacios et alJl2003|) . As 
these processes are not included in our models, we consider 
here two illustrative cases when prescribing Q(r): 

No mass loss is considered at any evolutionary stage, 
that is, the total angular momentum is assumed to be 
conserved. In this framework, the two assumptions for the 
transport of the angular momentum are: 1) either instan- 
taneous transport of angular momentum in the whole star 
(global conservation) which thus yields a uniform rota- 
tion, or, for sake of simplicity and illustrative purpose, 
2) local conservation of the angular momentum (shellular 
rotation) . Formally, the local conservation of angular mo- 
mentum between two instants (£2 > ti) can be expressed 
as: ^ 

-r?(ro) n x (m) = -r 2 2 {m) fi 2 (m) , (42) 

where m is the Lagrangian abscissa coordinate. Rotation 
in the convective core is considered to be rigid. As shown 
by Eq. gravity is modified by the radial term of the 
centrifugal acceleration A c (see Eq. l4*T|) . In Fig. ^ Ac is 
displayed as a function of the radius for both cases (uni- 
form and shellular rotation). The radial distance r = 0.1 R 
corresponds to a local maximum which is the limit of the 
convective core. This implies a modification of the local 
density, particularly in such regions. 

In Fig. [21 evolutionary tracks for a 1.8 M© model 
computed with three different assumptions: uniform rota- 
tion (Jg), shellular rotation (Jl) and absence of rotation 
(O = 0). Two evolutionary stages are labeled: the zero age 
main sequence (A, Aj and Al models), and on the main 
sequence (B, Bj and Bl models). Fundamental stellar pa- 
rameters of these models are given in Tabled During the 
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Fig. 4. First (a), second (b) order frequency corrections 
v = uj/2tt for 1=1 modes as a function of the radial 
order n. Finally, both contributions cji + UJ2 are displayed 
in panel (c). Empty and filled circles represent frequency 
corrections computed in the cases of uniform and shellular 
rotation respectively. 



main sequence evolution of the star, the rotation profile 
along the stellar radius mainly results from 1) the contrac- 
tion of the core, and 2) the expansion of the outer layers. 
Then the local moment of inertia of contracting shells of 
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Fig. 5. Weighted j/oi functions versus the radial distance r (normalized to the radius of the star R). Solid and dashed 
lines represent the / functions computed for the Bl model. Dash-dotted lines represent the / function (defined by 
Ea. I44JI computed for the Bq model. Dotted lines represent the rotation profile given by the radial function 770 (^) 
(the scale has been adapted for the sake of clarity). Panel (a) shows / computed for the n = 2,1 = 1 mixed mode. 
Respectively, panel (b) represents / computed for the n = 15, £ = 1 p mode. In both plots, the corresponding shellular 
rotation profile (same as Fig. \$ is also displayed. 



a given elementary mass decreases. Therefore, the local 
rotation frequency of those shells must increase in order 
to conserve the angular momentum locally. Consequently, 
the rotation frequency toward the surface decreases during 
the evolution of the star (see Fig. EJ. This figure shows a 
strong gradient of the rotation frequency toward the cen- 
ter: At a height of r = 0.2 R the rotation rate is fl — 1.5il s 
with f2 s the surface rotation rate. In our model, this means 
a rotation frequency around 14.97 yRz at r = 0.2 R, while 
a rotation frequency of 23.78 /J3z is obtained at the limit of 
the convective core (r ~ 0.1 R). Oscillation modes having 
amplitude in these regions (typically g and mixed modes) , 
are expected to be sensitive to shellular rotation. In prac- 
tice, observations provide the surface velocity (most often 
through the usim parameter). In order to reproduce such 
conditions, models Bl & Bq have been evolved so as to 
obtain the same rotation velocity at the surface. The first 
one (Bl) has been evolved considering the assumption of 
local conservation of the angular momentum. The second 
one (Bq) has been evolved under the assumption of uni- 
form rotation (global conservation of the angular momen- 
tum). 



5. Effects of shellular rotation on oscillation 
frequencies 

Oscillation frequencies of a rotating star are obtained from 
Eq. For the sake of clarity and shortness, hereafter, 
the following nomenclature will be used: luq and ojq for 
zeroth order frequencies, ui\ and for first-order fre- 
quency corrections and finally, u>2 and will be used 
for second-order frequency corrections. The L and G su- 
perscripts follow the same nomenclature adopted for the 
models Bl and Bq. Similarly, differences of frequency for 
a given mode (given n, I, m subscripts) between two mod- 
els will be denoted as Aujj = uj^ — uif for j = 0, 1, 2 associ- 
ated with zeroth-order frequencies, first- and second-order 
frequency corrections respectively. 

It is found that the behaviour of frequency correc- 
tions does not significantly change for low £ degree val- 
ues (£ < 3) when a uniform or a shellular rotation profile 
is assumed. We therefore focus on I — 1 modes in this 
section. 

In the following sections, the analysis and discussion 
of the effect of shellular rotation on oscillation frequencies 
include very high-order p modes. The reader should no- 
tice that such modes present a low probability of being 
observed (even from space) due to their expected small 
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Fig. 6. Differences between oscillation frequencies v = w/27r given in Fig. as a function of their radial order n. 



Such differences are displayed for the zeroth-order Aloq 



Alji = u{ - 
dotted lines. 



(represented by a continuous line); the first-order 



indicated by dashed lines, and finally for second-order corrections Aw2 = u>% — u>2 represented by 



amplitudes. Even so, from a theoretical point of view, in- 
vestigation of the frequencies of these modes - i.e. with a 
well defined (asymptotic) behavior - is found helpful for 
the global interpretation of the oscillation spectrum. 



5.1. Zeroth-order frequency corrections 

Perturbing and linearizing the hydrodynamical equations 
about a pseudo-rotating model yield an eigenvalue system: 
(see En. I A. 51 in Appendix El for more details). Our numer- 
ical oscillation code is adapted to solve this system, with 
the boundary conditions specified in Eq. IA.16tlA~19l and 
including a radial shellular rotation profile (Eq. • The 
resolution of this system yields the zeroth order frequency 
luq. The use of pseudo-rotating models described in Sect. 21 
provides zeroth order eigenfrequencies which include the 
effects of the spherically symmetric perturbation of pres- 
sure, density and gravitational potential. These perturba- 
tions are induced by the spherically symmetric component 
of the centrifugal acceleration (see SGD98). Eigenmodes 
of this system are still m degenerate. 



Oscillation frequencies are thus computed for the Bl 
and Bq models described in Sect. 21 In order to quantify 
the effect of shellular rotation on zeroth-order frequencies, 
A wo = oJq — lJq differences are calculated. These differ- 
ences are shown for £ = 1 modes in Fig. El (solid line) as a 
function of the radial order n of the mode. As can be seen, 
the effect of shellular rotation on zeroth-order frequencies 
is very small at low frequencies and increases with the 
radial order n. 

For high radial order p modes (n > 12), the differences 
Aug are larger than differences from higher-order correc- 
tions Awi,Acj2, (Sects. [3731 13.4JI and hence become the 
most sensitive to the effect of shellular rotation. 

Numerical tests show that the differences Aw mainly 
come from the differences in the pseudo rotating models, 
Bl 7 Bg; the contribution of the rotation profile (either 
uniform or shellular ) in the oscillation equations is negli- 
gible. 

As an order of magnitude, for ^ = 0, rt ~ 10— 15 modes, 
we consider the slope, v a [y = uj/2ir), of the variation of 
the differences (fo,ra+i,£=o,o — v a,n,l=0,o) with the radial or- 
der n (in the asymptotic approximation). We numerically 
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Fig. 7. (a) Same as Fig. but first- and second-order frequency corrections have been included. Symmetric solid line 
branches represent from top to bottom, the differences for m = — 1 and m = +1 mode frequencies respectively. For 
m = modes, differences are represented by a dotted line. Finally, panel (b) is equivalent to (a) but the the correction 
for near degeneracy is included. In both panels, the shaded region represents an indicative frontier between the region 
of g and gp modes (left side) and p modes (right side) . 



find that Av a = {v%-v%) ~ 0.175 /Hz for v a ~ 47/JIz i.e. 
a relative difference Av a jv a ~ 3.7 10~ 3 . This difference 
can be accounted for by the fact that our Bl and Bq mod- 
els have slightly different radii A R/R = (R L -R G )/R L = 
— 1.76 10~ 3 (see TableQJ. This induces a contribution to 
Av a which can be quantified as: 



2.6 ID-. 

Va 2 R 



(43) 



Hence in the present case, quite large differences, Aw , 
at high frequency, come from structural properties of the 
two models through differences in their structural vari- 
ables (pressure, density perturbations... ) rather than from 
differences in the wave propagation properties which ap- 
pear as additional extra terms in the oscillation equations 
(see Appendix El . 

5.2. First-order frequency corrections 

In Fig. Ok, the first-order frequency corrections u>i (Eq.EJ 
are represented as a function of the radial order n. Empty 



circles represent the frequency corrections obtained for 
the Bq model, i.e. in the case of uniform rotation. 
Respectively, filled circles are those obtained for the Bl 
model, that is, in the case of shellular rotation (i.e. assum- 
ing local conservation of the angular momentum) . 

Figure can be seen as composed of two domains: 
a first one (n < 5) containing g modes and mixed modes 
(hereafter gp modes region) , and a p -mode region (n > 5) . 
As expected, modes in gp region are more affected by 
shellular rotation than those in the p-mode region. 

In Fig. El the frequency correction differences, Au)\ = 
u>i — u>i , are represented as a function of their radial or- 
der n for t — 1 modes. These differences assess the effect 
of assuming a shellular rotation instead of a uniform ro- 
tation. In the region of gp modes, the effect of shellular 
rotation can reach up to 3 //Hz for a few modes. 

For high radial order n, differences AcJi, although 
smaller than for gp modes, are significant and remain 
quite stable around 1 fMz. 
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These results can be explained by analysing the con- 
tribution of shellular rotation to u>i corrections (Eqs.|HlE3 
fc lT2|l . which arises from: 

1) the effect of f2(r) on the zeroth-order eigenfunctions 
yoi and zq (corresponding respectively to the radial 
and horizontal displacements) appearing in Cl and Jo 
terms. Those eigenfunctions are obtained by solving 
the eigenvalue system taking into account a shellular 
rotation, fl(r), given in Annendix lAl fEns. IA.lt]A~4j) . 
Importance of this effect can be measured with the 
weighted function / as defined by 



f(r) 



4 2 
P T Vol 1 



(44) 



where yoi represents the radial displacement eigen- 
function. 

2) the explicit contribution of the rotation profile i]o in J 
(Eq. El • This latter term vanishes in case of uniform 
rotation. For high-frequency p modes, the horizontal 
component of the displacement is negligible (|yoi| >> 
00 ) and one has Jq ^< rjo > where < • > is a weighted 
average defined as: 



< ■ >- 



Jo P ri yli ■ dr 
lo P r4 Voi dr 



(45) 



In Fig. El the weighted function f(r) = pr^y^ and 
Vo(f)f(r) are then represented as function of the radial 
distance r for a high radial order p mode (Fig.^Ja) and for 
comparison in Fig. for a mixed mode. 

For gp modes, the eigenfunctions (as shown by / func- 
tions in Fig. QJO present an inner maximum near the core 
r/R ~ [0,0.2]. As the kinetic energy of these modes are 
large in core-close regions, large differences between the 
uniform and the shellular rotation cases i.e. in Awi are 
reasonably expected. 

For high radial order p modes, implicit effects of shel- 
lular rotation on eigenfunctions fi, are represented by 
dashed lines and can be compared with fa for a uniform 
rotation. Maximum amplitudes of such modes are located 
near the stellar surface, the eigenmodes do not differ much 
between the uniform and the shellular cases. The main 
contribution to the differences Awi = 1 pHz comes from 
the large increase of the rotation rate toward the center 
which is accounted for by 770 in J . Hence, negative values 
of Jo explain that wi absolute values obtained in the case 
of a shellular rotation are lower than those obtained when 
a uniform rotation is assumed as seen in Fig. QJi. 

5.3. Second-order frequency corrections 

Second-order frequency corrections u> 2 are found smaller 
than the first-order frequency corrections except at large 
frequency. They are presented in Fig. QJd for I — 1 modes 
as well as the sums (wi + lu 2 ) as a function of their radial 
order n. At low frequencies (n < 5), Aui 2 is comparable 
with the effect on second-order corrections (Aw ). 



In the p-modes region, investigations of the individual 
terms in Eq.EI reveal that X 2 and Y 2 are clearly dominant 
in that region. For given £ = 1 centroid modes, w 2 can be 
written as: 



uj 2 {£ = l,m = 0) = 



Q2 

U 



X, 



(46) 



In this expression, X 2 is found to be dominant mainly due 
to its wo dependence through the 3 C term (see Eq. [E 
Neglecting then Xf and X[, and considering X 2 ( 
l,m = 0), Eq. 021 reduces to: 



uj 2 (I = 1 , m - 



0)~-— 3c 
5 w 



(47) 



For m = ±1 modes, Y 2 behaves like X 2 but with opposite 
sign. Similarly to X 2 , the presence of Cwq in 3c makes 
Y 2 dominant face to Y\. Thus the shape of the m — ±1 
branches shown in Fig. QJd can be explained by the compe- 
tition between X 2 (positive) and Y 2 (negative), yielding: 



w 2 (£ = l,m^0)~-i— 3c 
5 w 



(48) 



More generally, for high radial order (£ ^ 0, m ^ 0) modes, 
one has (see also DG92): 



U! 2 • 



(A -3 m 2 ) n 2 



4A- 3 



w 



3c 



(49) 



Unlike the first-order frequency corrections, the second- 
order m ^ branches are not symmetric with respect to 
m = modes and the asymmetry for a given £ can be 
written as: 



^2 : m + U2,- 



2^2.0 = -- 



6m 2 n 2 
4A-3^o" 



3, 



(50) 



For high radial order p modes, 3c is given to a good ap- 
proximation by 



3 C ~ fJ < §2 > 



(51) 



where a 2 — ujq/(GM/R 3 ) < S 2 > is a weighted average 
of perturbations of structure (see Appendix El En. lB.12)i . 
Hence, 



to 2 (£= l,m = 0)~ - 



n 2 



5 GM/R 



■W < §2 > 



(52) 



and the increase of W2 with n clearly seen in Fig. QJ) arises 
from its wo dependence. In the -modes region, in con- 
trast to pure p modes, X\,Y\ (more specifically Xj , Yf) 
dominate over the structure rotational distortion effects 
X 2 ,Y 2 . In that case, both effects 1) and 2) mentioned 
above contribute in an intricate way. 

The combination of both first- and second-order fre- 
quency corrections in Fig. Q] is then found to be domi- 
nated by the u> 2 behavior at high frequency and as it was 
the case of first-order corrections, u> 2 present absolute val- 
ues higher than u) 2 . At low frequency (gp -modes region), 
a complicated behaviour results from both contributions 
(wi + w 2 ). 
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Fig. 8. Variation of Av — Auj/2tt with Aw = Awo+Aa->i + 
Aw 2 as a function of the rotation rate of the uniformly 
rotating convective core normalized to the surface rotation 
rate Q/fi s = 1.2. Effect of varying the rotation rate is 
shown for modes in gp- and p-modes regions. The plots 
focuse on i — 1, m = 0, ±1 modes for which Aw gets a 
maximum value (see Fig. uh)'- n = 3 modes in the gp 
domain (panel (a)), and n = 16 modes in the p-modes 
region (panel 6). 



Figure shows the frequency differences Au2 — lo^ — 
w G . In general A1V2 < AcJi over the whole range of ra- 
dial order considered. Contrary to the first-order results, 
the effect of shellular rotation on second-order frequency 
corrections (A1W2) is larger for p modes than for gp modes. 

For high radial order p modes, the second-order fre- 
quency differences arise from: 

CP 



A0J2 



GM/R 3 



A(tU < §2 >) 



with 



uj5 < § 2 > L 



< §2 > 



G 



(53) 



(54) 



A(w < §2 >) 

As shown in Fig. El in the p-mode region, Aloq increases 
with the radial order n. The effect of this increase is low- 
ered in Aw 2 by the change of < § 2 > which itself is mainly 
due to the derivatives of f2(r) (i.e. 62, dbzjdr). 



For gp modes, the difference in uiq between both mod- 
els is small (Fig.EJ and its contribution to Aw2 (Eq.44) is 
less important. In this case, second-order effects of cen- 
trifugal and Coriolis forces on the stellar structure are 
dominant through terms in # c proportional to: 



du 2 
dr 



ru 2 , 



rb-2 



(55) 



This is shown in equations describing the behaviour of 
the pressure and density distributions in presence of ra- 
dial shellular rotation (Eqs. 74-80 in DG92). Moreover, 
second-order effects of shellular rotation on the non- 
spherically symmetric component of the gravitational po- 
tential $22 are found to be negligible and 112 ~ (1 +772) in 
Ea. IB. 81 Therefore, Au>2 differences are dominated by the 
effect of ft(r) (ryo) and of its derivatives (62,...). 

As it was the case of first-order corrections, ojt? present 
absolute values higher than u>2- 

5.4. Full non-degenerate frequencies 

Finally, the total effect of differential rotation up to second 
order on the frequency Auj is displayed in Fig.EJi. 

The largest effects of shellular rotation on gp modes 
(up to ~ 3.4 /Hz) arise for m = ±1 modes and are dom- 
inated by Aoj\. In this region, the effects on zeroth order 
frequencies are much smaller. In contrast, for the p-modes 
region Auj are dominated by Awo increasing with the ra- 
dial order n (see Sect. 15. 1|) . Numerically, we find differences 
up to 0.5 pHz for m = +1 modes. For m = — 1 modes, ef- 
fects up to 3 pHz are seen to be of about the same order 
than those found for m = modes. 

As already mentioned, the efficiency of the transport 
of angular momentum is modified by rotational mixing 
during the evolution of the star. The effect of shellular ro- 
tation on oscillation frequencies is expected to be reduced 
(in preparation), due to the smoothing of the variation 
of Q/Q 8 at the edge of the convective core (see Fig. EJ. 
In order to simulate such a situation, theoretical oscilla- 
tion spectra for the model Bl have been computed with 
core rotation rates in the range of fl/tt s = [1.5, 1.2]. As 
expected, the major influence of such variations is found 
for gp modes. The lower the rotation rate in the core, 
the lower the effect of shellular rotation in this region 
(see Fig. |Ht)- For the minimum rotation rate considered, 
\Au>\ ~ 1.4 pHz for m = — 1 modes, and |Aw| ~ 1.5 /Mz 
for to = +1 modes is found. For p modes (see Fig. Et), 
differences remain close to those obtained for fi/f2 s = 2.5 
(used along the present work). 

5.5. Correcting for near degeneracy 

Near degeneracy is taken into account as described in 
Sect. The oscillation frequency of near-degenerate 
modes are computed according to Eq. EHI & 1221 or 
Eq. d & E3 In Fig. 03, oj l and uj G are represented for 
(i a ,h) = (1)3) degenerate pairs. Near degeneracy is con- 
sidered here for modes with \a a -a h \ < fi s / (GM/R 3 ) 1 / 2 - 
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Fig. 9. Effect of the second-order near-degenerate coefficients C ab /uj a and C ab /ujb as a function of the frequency 
(y = uj/2tt) for degenerate pairs (£i,£ 2 ) = (0,2) (left panel) and (£i,£ 2 ) = (1,3) (right panel). Empty and filled circles 
represent degenerate coefficients in the cases of uniform and shellular rotation respectively. 



0.25 (i.e.\v a — Vb\ ~ 10 yRzhere) corresponding to the 
surface rotation frequency of Bq and Bl models (see 
Tabled}. As can be noticed, for both models, modes in the 
gp region clearly are the most affected by near degener- 
acy, with a few m = ±1 near-degenerate modes within the 
range of n < 8. This result was expected since such modes 
are likely closer to each other than p modes. In order to 
quantify how the oscillation frequencies are modified by 
near degeneracy, we define the coefficients and C^? 
composed for first- and second-order contributions, where 
ab subscripts represent the two modes affected by near de- 
generacy. According to the definitions Ecl UM!?! " + " sign 
is for mode a and "-" sign is for mode b. In the case of 
near-degenerate modes with same degree £, the first- and 
second-order coefficients are respectively defined as 

C ab ' = wi - wi ± C K ab ' = u 2 -u 2 -—. (56) 

In the case of two degenerate modes with different degrees 
I a = £b±2, the first-order coefficient C^' = (see Sect.EJ 
and 

C ab =U2-U 2 ± — (57) 

The full expression for the adiabatic oscillation frequencies 
can be written as: 

w = f "0 + [«i + C#] + [lo 2 + C$] if 4 = lb 

\lU + LO! + [U 2 + C^] if4=4±2 



In Fig. for degenerate pairs (£ a ,£b) — (0,2) (left 
panel) and for (£ a ,£b) = (1,3) (right panel), the effect of 
near degeneracy is assessed by C^ b /u>o as a function of the 
frequency for the cases of uniform and shellular rotation. 
As commented above, for such degenerate pairs, only the 
second-order degenerate coefficient is nonnull. It can be 
shown that the dominant terms correspond to w 2ah and 
^2 'ab fEas. KTTH and IC.llj) . The first one takes into ac- 
count the distortion of the structure of the star through 
the structure quantities d\ and its derivative d 2 (Eq.C20). 
The second one, wJab corresponds to the toroidal compo- 
nent of eigenfunctions. For the range of frequencies treated 
here, the global effect remains small for high-frequency p 
modes with respect to what is found in the -modes re- 
gion, where w^ab an£ l ^Jab are °f t ne same order. This 
explains that differences between differential and uniform 
rotation are mainly found for gp modes. Such variations 
can represent a few percent of the value of frequencies not 
corrected by near degeneracy. 

For m = modes, compensations of dominant terms 
turn out to yield a smaller effect of near degeneracy for 
m = modes seen in Fig.|Hl(left panel), compared to those 
with m 7^ (right panel). 

For modes in the high frequency p-mode region, the 
near degeneracy effects shown in Fig. can be studied 
in the asymptotic approximation by analyzing individual 
terms of Eq. IC.14I in Appendix IC.3I This analysis reveals 
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that for large frequency p modes the small separation cor- 
responding to v n -i,i — v n .i+2 is small enough with respect 
to the rotation rate of the star to degenerate all (£, 1 + 2) 
pairs. In this context (asymptotic regime), near degen- 
eracy contributes with a slowly-increasing term with ra- 
dial order, dominated by wi? ab (see details in Sect, IC.3~T)l . 
which is given to a good approximation by 

Cl 2 24 

UJ < i>2 > , 



W 2,ab 



(58) 



LUq < § 2 > . 



(59) 



GM/R 3 7V5 
for (£ a , £b)=(0,2) degenerate pairs and 

D n 2 12 J2 

W2 ' ab ~ GM/R 3 TV 7 

for (£ a ,£b)=(l,3) degenerate pairs. For high-frequency p 

(2) (2) 

modes, C ab increases as Qq hence the ratio C ab ' /u>o re- 
mains nearly constant. 

Fig. also reveals that the behavior of the cj^ coef- 
ficients differ when differential or uniform rotation is con- 
sidered. These differences are found significant for m =^ 
degenerate modes and are mainly caused by differences at 
zeroth order (uj in Eos.l58lfc l59)l . 

For modes with A£ — (for which C^jJ ^ 0), further 
information can be obtained through the second-to-first 
order coefficient ratio 



9(2,1) 



c 



(2) 
ab 



c. 



(i) 



(60) 



In the present case, such degenerate pairs (which repre- 
sent only a few marginal cases) are found for 1 — 2 and 



t = 3, for which cj^ ~ C^j . When considering differen- 
tial rotation, mixed modes (n = —6. .2) present second- 
order degenerate coefficients slightly higher (from 1% to 
20%, in absolute value) than first order ones. Although 
the correction for near degeneracy is thus less important 
in the case of shellular rotation, it cannot be neglected. 



6. Discusion and Conclusions 

In the present work, adiabatic oscillation frequencies of 
stellar models with intermediate mass have been com- 
puted, which include effects of rotation. Following the per- 
turbation formalism of DG92 and SGD98, we have built 
a numerical code which provides adiabatic eigenfrequen- 
cies corrected up to second order in the rotation rate f2 
allowing a radial dependence of ft = fi(r). Two particular 
cases for the transport of the angular momentum in stellar 
interiors are investigated. On the one hand, an instanta- 
neous transport of the angular momentum thus leading to 
a uniform rotation profile. On the other hand, a shellular 
rotation profile f2(r) is assumed and derived from the as- 
sumption of local conservation of angular momentum. Our 
study has been focused on 1.80M Q equilibrium models of 
typical 5 Scuti stars with a surface rotation of 100 kms -1 . 

As mentioned in introdution, seismology of 
A-F rapidly rotating stars faces many difficul- 
ties, among which, the most severe probably 



is m o de identification (s ee reviews Gou pil et alJ 
120051 : iDupret et"ail I2005L for S Scuti star and 7 
Dor stars respectively) We stress that the present 
study does not intend to solve the problem of 
mode identification which starts now to be han- 
dled with confidence using various techniques (see 
Handler] I2005L and references therein). In par- 
ticular, diagnostics based on multicolor photom- 
etry l[Ciarrido et al1ll990 ^ are nowadays employed 
for mode identification. In the framework of the 
forthcoming space mission Corot, the colour pho- 
tometric data from the exo-planet camera will pro- 
vide very useful information for this task. In any 
case, whether mode identification is successfully 
performed or not, calculation of models and asso- 
ciated oscillation frequencies as close to reality as 
possible is necessary for a one to one frequency 
comparison between observed and computed fre- 
quencies (i.e. spectrum of frequency differences or 
frequency histograms). At that level, for establish- 
ing possible seismic diagnostics, it is important 
to be able to compute frequencies with different 
physical assumptions and assess their respective 
effects. The present study has been made along 
this viewline. 

To estimate the numerical and theoretical uncertain- 
ties of the computed eigenfrequencies, as wel as the va- 
lidity of the adopted perturbation technique is not an 
easy task. Here we have performed a one-to-one eigen- 
frequency comparison between eigenfrequencies computed 
with the present numerical code and those computed with 
the Warsaw oscillation code (restricted to second-order 
perturbation frequency corrections). This excersise also 
implies the comparison of the different perturbation tech- 
niques employed in each code. The input model for both 
codes is the same, Bg, with uniform rotation (Tab. QJ. 
The results show that differences less 0.15 [Mz are most 
oftenly obtained and never larger than 0.5 [Mz. Both oscil- 
lation codes are different in their numerical schema, treat- 
ments of boundary conditions and other technical com- 
putation characteristics. The adopted perturbation tech- 
nique is also different as, in one case, the first-order fre- 
quency correction is explicitly calculated (present code) 
whereas in the other case, the first-order effect of rota- 
tion is implicitly included in the eigensystem. The results 
of the comparison hence give some estimate of the preci- 
sion of the frequencies computed here and give us some 
confidence in their calculation. 

Note also that triple-mode interaction occurs between 
three modes such as £ a — £b — £ c — 2 or £ a — £ — b — 
2 = £ c — 2 in the adopted perturbation formalism here, 
although only occasionally; these triple-mode interactions 
correspond to a double-mode interaction in the SGD98's 
formalism since in this last one, the case l a = lb near 
degeneracy is implicitly included. 

Important effects - up to 3 /iHz - of shellular rotation 
compared with a uniform rotation are observed in the low 
frequency region where gp modes are encountered. This is 
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expected as the rotation profile rapidly varies in the inner 
layers at the edge of the convective core where gp modes 
have large amplitudes. 

In the case of higher frequency p modes, frequency dif- 
ferences up to 3 [Mz are predicted between a uniform and 
a shellular rotation. This is due to the structure deforma- 
tion caused by the centrifugal force, which mainly affects 
the zeroth-order oscillation frequencies. 

It is found that, the impact of reducing artificially the 
core rotation rates of models by 50%, principally affects 
gp modes for which the effects shellular rotation are dras- 
tically reduced to 0.6-1 /xHz. In addition, the frequency 
changes due to a shellular rotation compared to a uniform 
rotation are found much larger for these intermediate ro- 
tators than third order (ft 3 ) effects (around 0.04 /J3z for 
the region of g and gp modes, and around 0.02 /J3z for 
p modes, for a 100 krns" 1 star model). This represents 
around 3.3% and 2% of effects found for the core-reduced 
rotation models commented here-above. 

Note evenso that such effects can even be 
larger for different kind of stars (Suarez et al., 
in preparation). Observational evidences for this 
have already been obtained for the early-B type 
stars HP 129929 and HP 29248 l|Aerts et alJ l2003: 
IPamvatnvkh et al.l l2004. respectively). 

For the models and rotation rates investigated here, 
quite a number of modes are near degenerate and a double- 
mode interaction had to be computed for these modes; on 
the other hand, we find that only a few modes require a 
triple-mode interaction i = 0, 2, 4 or I = 1, 3, 5: one set at 
high frequency n = 9, 10, 11 and two or three sets at low 
frequency (in the fundamental mode vicinity). For these 
triple near-degenerate modes, we encounter all possible 
cases of the two closest modes among the three being ei- 
ther ^ = 0,^ = 2or^ = 0,^ = 4or^ = 2,4; this is the 
same for i = 1,3,5 modes. The effect of taking into ac- 
count triple interaction instead of double interaction can 
be quite different depending on the configuration of the 
3 modes (closeness of the frequencies). We find that the 
effect can be occasionally amount up to a few /Hz. 

In view of our different tests as indicated above, we 
conclude that if one looks for seismic signatures of shel- 
lular rotatio n smaller than 0.5 pft z, we have to resort to 
SGD98 and iKarami et alJ ll2005h formalism and include 
triple-mode interaction whenever it is necessary. 

For higher rotation rates, perturbation techniques are 
no longer useful, if even still valid, when the effects 
of rotation are large enough that triple-mode interac- 
tion must be systematically taken into account for most 
modes; this situation indeed would suggest that higher 
than triple-mode interaction ought to be included. Then 
no n-perturbative techniques su ch as the one developed 
by iLignieres fc Rieutordl J2004T) and Rees et al. 2005 (in 
prep.), for instance, will have to be used. 

In the framework of the forthcoming sp ace mission 
COROT jBaglin fc The COROT Teamlll99^ . precisions 
of 0.1 [Mz on mode detection, and around 0.5 [Hz on split- 
ting resolutions are expected to be observed. With such 



an accuracy, effects of shellular rotation are likely to de- 
tectable, provided numerical eigenfrequencies reach this 
level of precision. Finally, we stress that it is possible to 
remove contaminations of high-order (second and higher) 
effects of rotation from the rotational splitting in order to 
recover the true rotation profile w ith usual inversion tech- 
niques: this has been out lined by jPziembowski & Gounill 
11994 iGounil et alJl2004T) and will be described in more 
detailed in a forthcoming paper. At low rotation rates, 
this task would become even easier, since, in this 
case, the effect of differential rotation on oscilla- 
tion frequencies is dominated by structure terms 
in zeroth and first-order frequency corrections. 
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Appendix A: Oscillation frequencies of a 
pseudo-rotating model 

Considering an effective gravity, g e ff fEn. 139(1 . and a rotation 
profile, £l(r) (Eq.0, the following dimensionless quantities are 
defined as in SGD98: 



2/oi 



1 i j P'\ 

2/02 = {</>'+ —) 

Seff r p' 



2/03 = 
2/04 = 



4>i 



(A.l) 
(A.2) 

(A.3) 
(A.l) 



9cffr 

1 d<j)l 
g e a dr 

Following lUnno et alJ fl9891. we obtain a linearized eigenvalue 
system: 

dyoi . „ A 
x—j— = A - 3yoi + ~ — j V° 2 
dx GrO-g 



= (C>f ' " A * )V01 + {A * + 1 ~ Ux)V02 ~ A * V0S 

X ~dx~ = ^ C/ x)2/03 +2/04 



(A.5) 



"2/04 r , , .1 

X ~dx~ = 1 - a I V(n + 9< - y ° 2 _ yo3 '\ + yo3 _ U * yoi ' 

with x — r/R and R t he stellar radius, where the following 
classical notations from llJnno et all l|l989h are used 



^ 1 dlnp dlnp 
r d In r dlnr 



V = 



U = 



dlnp 
dlnr ' 
dlnM r 



v a = l 



din r 

A = £(£+1), 

and with U y defined as follows: 



(A.6) 

(A.7) 

(A.8) 
(A.9) 



U x = U + x- (A.10) 
We also define, as in SGD98, the following variables 

c=(;) 3 ^ a = T^- (a. ii) 

2 i^o 2 -R 3 A c 



(A.12) 

(A.13) 
(A.14) 

where A c is the radial component of the centrifugal accelera- 
tion defined in Ea l41l M, m r respectively are the stellar mass 
and the mass enclosed in the sphere of radius r. This system 
is solved with the following boundary conditions: 



G M 

ffeff V dr 
A = V g (yoi - yo2 +2/03) , 



yo2 + 2/01 4- = 0, 3y i + yo4 = (£ = 0) 



(A.15) 



2/01-^02^=0, yo4-%)3 = (i^0) (A.16) 
at the center of the star and, 

2/01 = 1 (A.17) 
2/04 + (£ + 1)2/03 = (A.18) 
A 4 + Col \ 



2/oi ^1 + 



A 



VCa'i V ) y ° 2(1 VCol> + 



J/03 ( 1 + — ) = I) 



(A.19) 
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at the stellar surface. The resulting eigenvalues correspond to y T — 1l f (1 + ^q) 2 

the oscillation frequencies of a pseudo-rotating model (with Io J a 

their associated eigenfunctions). „ 

+ (yoi + (t + l)z ) 2 



(yoi - £z Q ) 
- 1 1 



,£ + 2 



1 



I U 2 



+ 1 4(^ + l) 2 
po r 4 dr 



(B.4) 



A.l. First-order perturbed eigenfuctions 

2 f R r 

Considering dimensionless variables equivalent to Ea. IA.lT4A. 41 Y 1 = — (1 + 770) j^oyi + 

with first-order perturbed quantities (£i,r, 4>'i, Pi), and the ' 11 

zeroth-order solutions obtained from Eq. IA.5I first-order per- 
turbed eigenfuctions can be calculated from the following sys- 
tem: 



Z1V01 ~ VoiVi 



(A - l)zi2 par dr 



dyi . „ A 

x— = Ai - 3yi + — — j y 2 
dx C r o"o 

+ (yoi + z )(l + 770) - (770 + (Ji)Az 

x ii = {Cra ° _ A * )yi + {A * + 1 " Ux)m ~ A * m 

+ (ffi + r]o)yoi - (1 + Vo)za 

^ = ( 1 ~ t/ x)y3 + y4 



(A.20) 



i(<7i — l) 2 . 



(B.5) 



(B.6) 



Furthermore, X 2 et Y 2 contain the effect of deformation of 
the star due to the non-spherical component of the centrifugal 
acceleration. They are given by: 



X 2 = 



A 



(A.21) 



4A-3 



Y 2 



4A - 3 



(B.7) 



V 



[AV + V 9 y 2 - V ff y 3 ] + Ay 3 - U x y t 



dx 1 — a 
where Ai = V g (yi - y 2 + ys) and 

yo2 

zo = Crf 



with 



2Io Jo 



(A.22) 



The horizontal component of ^ 1 can be written as follows: 



2/2 



l+?7 o, 
A 



-yoi 



l+7?0 

A 



01 20, 



(A.23) 



where 01 = Cl — Jo represents the first-order correction of the 
corresponding eigenfrequency. 

Appendix B: X{ and Yi expressions in Saio's 
notation 

The terms Xi and Yi are constructed from lu 2 ,lo 2 , ui 2 and, cu^ 
(see Sect.EJ. It is convenient to split X\ and Y\ as follows: 

Xi = x[+x? 

Yi = F/ +Y 1 r + Y 1 P + Y? 



[Col [(2yoi 2 + 3*g)(r^ + (4 - t/> 2 ) 

+ 2zqAu 2 + zoyoi(*i + 6m 2 )] 

+ yoi (A - yoi (A* + V g ))Uu 2 — yoiyo3(*i + 2A7J2) 

- yoiyo4(3r^+ U2 (10-3C/)) 

- 2z yo4A7i 2 - [yoi 2 (Ccro + 3 - U) + ^ 2 C7a 2 (A - 3) 
+ 3z yoi - 2y iyo4 - 2z yo3(A - 3)] b 2 

2 db 2 1 4 

- yoi r -^rj Po r dr, 

with it2, 62, »?2 defined as in DG92: 

112 = r t§ + 5 (1 + 7?2) 
, 1 dr) 2 
62 = 3^ 

?72 = (1 + 7?o) 2 - 1 • 

B.l. High radial order p modes 



(B. 



(B.9) 



These two terms include both, the toroidal component of the For high radial order p modes, the horizontal component of the 
first order perturbation of the displacement eigenfunction and fluid displacement is much smaller than the vertical one and 



an effect of inertia. The different components of Xi are given I yoi I >> \zo\- The modes are concentrated toward the sur face 

and yo3,yo4 ~ 0, U ~ 0, C ~ 1 , A ~ V g y i. Then Eq. |Bj 
becomes: 



by the following analytic expressions: 

X ' = /0 A 1 r -\r I fe 2[2(A - l)yoi 2 - Ayoizol Por 4 dr 
' - 6)Io Jo L J 



(8A- 



X\ — — . 

1 o Jo 



(B.l) 



+ ^o) 2 [(yoi - izo) 



a (£+l)(£+2) 
4(^ + l) 2 - 1 



Wo I y?n{c^mr^ + ^ 2 )- b2 ] 
- (Ca 2 +3)b 2 -r^}p r 4 dr. 
For large (To, one further has 



(B.10) 



+ (2/01 + (£ + l)z ) 



2 £(£- 



u 2 - 1 

Respectively, the components of Yl are expressed as 
1 



]por 4 dr. (B.2) j c ^ a 2 l. j* y 2 n { (r ^ + ^ _ ^ j ^ r 4 dr (fi n) 



1 

2/0 7 



?72 [yoi + (A - 2)z - 4y i^o] Po r dr (B.3) 



As 



then 



2 4 , 

y i Po r dr, 



~ O"o < §2 >, 



(B.12) 



where <> is a weighted average defined in Eq. Hence 

J c (thereby X2,Yi) increases with a 2 i.e. with increasing fre- yoi,aVib + AaZoi.a^ib 

quency. 

B.2. Radial modes 

For radial modes, £ = 0, one has X2 = Eq. 03 simplifies as: 
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with C.2. Detailed expression for u^lL an< ^ ^2 ab 

/ du2 , \ . „_> In the case of cJoTk, its components are expressed as: 

< S2 > = < [r—j — h 4u 2 - 6 2 J >, (B.13) • 

^ 2 P ab = / Vpo dr (1 + 770) 2 { (C.8) 

Ja^O Jo L 

H.aJ/lb + A a Z01,a^lb 

— (yoi.a^lb + Z01,aJ/lb + 201,a«lb)| , (C.9) 

and 

W^ab = / tA P0 dr(l + f?o) 2 { 

^2 /•« / 4 s ^ a A (n+iTb+iA^+i +fi_ifb_iA J e fc _i) (CIO) 

^2 = / ?/oi 2 ( 62 + -(1 + r?o ) ) po r 4 dr (B.14) - / . \ x , . \\ 

7oUoio V 3 / +°tah+2 [n+in+l Ab+lJ + <5£' a ^_ 2 (Tt-iTb-l A(,_lJ j, 

where t},+ i and t&_i are defined as follows: 

Appendix C: First- and second-order corrections & t \ tn-ti\ 

c , ^ , . n+i = - , . (yoi.b -4z ,i,)7ft+i (C.ll) 

for near-degenerate frequencies (.*& + ^ 

The different terms which contribute to the first- and second- ^-i = 7" (jfoi.b - (4 + l)zo,&) 76 ■ (C.12) 
order frequency corrections in presence of degeneracy are col- 

lected in the N and ^ given by Eq. EU These contributions Finall y< the mertlal term of the correction is given by 

are : 1 2 ^ 2 f R 4 2 f 

^2,ab = -m or _ / r p dr(l + r/ ) ^ 

I iJaWO Jo 

-j-(£o,aW£o,b) (CI) (yoi, a yoi,b + A b+1 z 0ia z ,b) 

~ [«0,»l K I tf.b) + (€ ,al 3C I ^b>] ~ 2 (»>1.»«0,» + lft>l,b*0,a + «0,.*0,b) } • (C.13) 

W 2 ,ab+wjab (C2) 



Wl.ab 



PT 
< jJ 2,ab 



^2,ab = 57T-p <^0,a|— (^2 - f>2 (J ) I € .b) (C3) 

^2,ab = 2 - o j [(io,J m2 ^ 2 ~ 2m ^\io,b)\^ (C4) torized as follows 

where 7 a fEa . II IB is the normalization term for the mode a, 

la = « ,.l€0,.>, (C5) 



For the correction u;3? a b, the angular dependence can be fac- 



k>2,ab = ®2ab Zab, (C14) 

WO 



where Cbat and Z a (, represent the angular and radial compo- 
nents respectively. Selection rules then exist and are imposed 
by: 

where the scalar product definition is defined as in Eq. (12) ^ f y*y pi 9)&Q, (C 15) 



of DG92. £ >1 k and £ x k represent respectively the poloidal and , 

toroidal components of first-order eigenfunction correction (see where ^ = sin9ded ^~ is the elementary solid angle and 

Eq 26 in SGD98 for explicit expressions). The operator DC is ^ ig thfl spherical harmonics y Using the sec0 nd-order 
defined as: 



3C = mfl - Qfle z x . (C.6) 



Legendre polynomial: 



3 2 1 

Definitions of the operators &2,p2 can be found in DG92 and 2 2 2' ^ 

SGD98 - we obtain 

_ 3 \ , 2 .2 1 

C.l. Detailed expression foruJi^ 



Q2ab = - [(7a + l +ll \- ^)3ab 

+ / ya'7b+lSab+2 + 7i>7a + l<5af)-2 I , (C17) 



The use of the definitions of £2 and 3C on Eos. lC.ll4C.3l vields 
the following expanded expressions: with 

mfi r [*_ 4 J „ s (■ 7 b =\/^ and y f = 1 ~ '" . (C.18) 

Wl.ab = 



nfl t / % , „ ,f 7fc = V ^6 and J £ = — -. 

~j-$l a l b r 4 p dr(l+r ?0 ) U 2 -l 

la Jo *• 

1 . \ On the other hand, the radial component Z a b can be writ- 

-(yoi,ayoi,b + A b zo, a zo,b) ten ag SGD98: 

+ (j/01,aZ0,6 + Z0,aV01,b + Z0.aZ ,b) k (C7) l f R ( 1 

Zah = Wl V P° Ar { d i E i+ d2E2 + r ^E 3 +r 2 b 3 E4,j, 
where &i a i b is the Kronecker symbol. (C.19) 
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where where 

a „ dd! , 1 2 d 2 V2 ,„ , Aj-lM + l), (C.32) 
d2 = r lF 63 = 3 r ^> < C - 20 > with i representing the modes a and b. 
with 62 defined in B9. The non-degenerate expression is a par- 
ticular case when a = b and it can be found in DG92. /-oiiri j-i j j 

C.3.1. High radial order p modes 
Ihe explicit expressions of ha, E%, E3 and £/4 are given 

below. For the sake of shortness, the following nomenclature is Following the same approximations for modes concentrated to- 
used in DG92: wards the surface as in Sec. IB.ll and using the following ap- 



_ „ „, _ „ tit _ „ tr< m\ nroximations ffrom the definitions Ens. IC.22l4C.24ll 

y = yoi v = y 3 W = y 4 , (C.21) ^ v ' 

and the following short expressions have been defined: 1i ~ 3/j (C°"o,j + 4) , Aj ~ V g yj, Sj ~ yj j = a,b (C.33) 

Qj = 2/j (Co-o.j + 4 ~ U) - A^o.j - Wj (C.22) En. ETT91 becomes for large cr : 

= ^(yi-iftKw+tfj) (C.23) R ^ 



pTiAj Zab ~ — Cao^4:U2 - b 2 jy a y b p r 4 dr , (C.34) 



.-. - . (C.24) /a 

rpQcft 

f or j — a yy e then have- where we have used definitions given in En. IB. Al and En. lC.2f)l 

This expression is the same as the one in Bll when a = b. 
Using the definition of weighted average given in Eq. 0^1 
Ek = E k , a b + Ekfi a k = 1, 4 (C.25) with y a yb instead of j/01, Eo. lC.34l can be rewritten as : 

Z a b ~ CTq < §2 >ab ■ (C.35) 

3 11 

Ei a b = --^qaVb + ^VaVb^i -q a Xb - ^ybS a if) D 

2 2 2 tmally, w 2 ab (Eq.C14) for high radial order p modes can be 

+ (Cffo,„ + 4)y a A 6 - A 6 Wa written as 

+ I U - - ) - A„v jfc w D ab _ — — Qaab w < S2 > ab ■ (C.36) 



A a z a . a W b + 3q a z , b + -zo.a^ct^ 





- A a z 0>a y b f | - Co-o, a + f/ -4j - 


■ A a ZQ,aX b 






- ZO.ay&S^f, +ZO,a2/Q2,b(A+ 6 — 3)(£7 


-4) 










(C.26) 


E2ab 


1 2 1 


- 2)y a s b 






1 , 2 

- VbW a + -(A ab - 3) Ccr 0ta Zo,aZo,b + A b z ob y a 






- (Aaf, - 3)(z ,by02,a) ~ ^S a A b Z 0b , 




(C.27) 


Ez a b 


= -2 C<T 0,a(Vayb + (Atb-3t)Z0,aZ0,b) 






~ ^Vo.s b {U + 6) - 32/f,2o,a — 2"(^- a " 


- A b )ybZO,a 






- (3 - A+ b )zo, a v b + y a W b 








1 fdlnT 1 \ 




(C.28) 


Ei = 


1/ 

= — 2^ aSb + yi" 8 ")' 




(C.29) 


with 










= 6 + J7({7 - 2(A* + l/ g ) - 3) + (3 - 


U)Co-l >a 




J 2a 


= (A+ -3)(L/-4) + A+(C7aL + 4 


-U) 






= (A+ -3)C7a 2 , a (6-t/)+6A+ 




(C.30) 




= (l-U)(4-U)+6 + U(U-3) 






A+ ft 


= I(A a+ A 6 ) 






A 6 ~ 


= - (A* - A a ) = -A~ b , 




(C.31) 



